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Abstract. A few years ago we identified a recursion that works di¬ 
rectly with the gaps among the generators in each stage of Eratosthenes 
sieve. This recursion provides explicit enumerations of sequences of gaps 
among the generators, which sequences are known as constellations. 

Over the last year we identified a discrete linear system that exactly 
models the population of any gap across all stages of the sieve. In August 
2014 we summarized our results from analyzing this discrete model on 
populations of single gaps. This paper extends the discrete system to 
model the populations of constellations of gaps. 

The most remarkable result is a strong Polignac result on arith¬ 
metic progressions. We had previously established that the equivalent 
of Polignac’s conjecture holds for Eratosthenes sieve - that every even 
number arises as a gap in the sieve, and its population converges toward 
the ratio implied by Hardy and Littlewood’s Conjecture B. Extending 
that work to constellations, we here establish that for any even gap g, if 
p is the maximum prime such that p^\g and P is the next prime larger 
than p, then for every 2 < ji < P — 1, the constellation g, g ,..., g of 
length j i arises in Eratosthenes sieve. This constellation corresponds to 
an arithmetic progression of ji + 1 consecutive candidate primes. 


1. Introduction 

This paper is a technical note that builds on our previous paper, Er¬ 
atosthenes sieve and the gaps between primes [8|. As such, we will assume 
familiarity with the notations, models, methods, and results of that work. 

A constellation among primes 03] is a sequence of consecutive gaps be¬ 
tween prime numbers. Let s = aid 2 • • • a*, be a sequence of k numbers. Then 
s is a constellation among primes if there exists a sequence of k + 1 consec¬ 
utive prime numbers PiPi+i ■ ■ ■ Pi+k such that for each j = 1,... ,k, we have 
the gap p i+j -Pi+j- 1 = a 3 . 

We do not study the gaps between primes directly. Instead, we study 
the cycle of gaps Q(p #) at each stage of Eratosthenes sieve. Here, p # is 
the primorial of p, which is the product of all primes from 2 up to and 
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including p. G(p#) is the cycle of gaps among the generators of Z rnodp^. 
These generators and their images through the counting numbers are the 
candidate primes after Eratosthenes sieve has run through the stages from 
2 to p. All of the remaining primes are among these candidates, but many 
of the candidates do not survive further stages of the sieve to be confirmed 
as primes. 

At times we will speak of constellations corresponding to gaps between 
primes. We do this for intuitive convenience, to motivate why the particular 
constellations we study may be of interest. The results here are for popula¬ 
tions within Eratosthenes sieve, and the constellations correspond to primes 
only to the extent that they survive the sieve. 

In pj] we identified a discrete dynamic system that exactly models the 
populations of gaps from the cycle of gaps G(pk#) in one stage of Eratos¬ 
thenes sieve to the next, G{pk+ 1 # )- We use the eigenstructure of this system 
to obtain asymptotic ratios between the populations of various gaps in the 
sieve. 

With the following theorem, we generalize the discrete dynamic system 
for gaps to model the populations of constellations across stages of Eratos¬ 
thenes sieve. 

Theorem 1.1. Let s be a constellation of length j\ and of sum |s|. Let po 
be a prime such that |s| < 2p\. 

Let nj(s,p #) be the population of driving terms of length j for s in G{p*), 
with j\ < j < J and p > po- Then 

(1) n s (p k *) = M ji:J \ pk ■ n s (p fe _i # ) 

= R ■ ^ji:j\ Pk ' L ■ n s {p k -i*). 

in which L is the upper triangular Pascal matrix, R is the upper triangular 
Pascal matrix with alternating signs (—1)* + T and A is the diagonal matrix 
of eigenvalues (p k — j — 1) for j\ < j < J. 

As with our work on gaps, this explicit enumeration has two limitations. 
The first limitation is that all of the populations grow super exponentially 
by factors of (p — j i — 1). To address this, we normalize the dynamic system 
by dividing by factors of (p — j\ — 1), which makes the dominant eigenvalue 
equal to 1. 

The second limitation of Theorem 1 1.1 1 is the condition |s| < 2p\. Because 
of this condition, to exactly model the population of a constellation we have 
to calculate its initial conditions in G(p#), for p very close to |s|. But G(p#) 
consists of <f>(p#) gaps, which becomes unwieldy around p = 37 : 
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p 

length of G(p*) 

31 

3.06561 MO 

37 

1.10362 E12 

41 

4.41448 £13 

43 

1.85408 £15 


The eigenstructure of the dynamic system in Equation (H|) enables us to 
extract the asymptotic behavior of the population for constellations, if we 
can derive the initial conditions for a cycle of gaps Q{p#) under much weaker 
conditions than those of Theorem 11.11 

Extending the analogue of Polignac’s conjecture for Eratosthenes sieve 
[8], we are able to establish a strong Polignac result for arithmetic progres¬ 
sions: that for any g = 2 n, every feasible repetition s = g,... ,g occurs in 
Eratosthenes sieve. To the extent that these repetitions survive the sieve, 
they correspond to consecutive primes in arithmetic progression with gap g. 


2. Modeling the populations of constellations 


2.1. Summarizing the model for gaps. In [7] we developed a dynamic 
system that exactly models the populations of gaps through the stages of 
Eratosthenes sieve. We restate the model briefly here. 

Let g be a gap (a constellation of length 1) with driving terms up to 
length J. Let po be a prime such that g < 2p\. Then 

n g (Pk*) = M\ pk ■ n g {p k -i*) 

0 ••• 0 

2 0 


Pk ~ 2 1 

0 p k - 3 


J- 1 

p k - J -l 


■ n g (p k - 1 *) 


with 


= R ■ A ip k ) ■ L • n g (p k - 1 # ) 

A (pk) = diag [p k - 2, p k — 3,... ,p k — J — 1] . 


Since the eigenvectors do not depend on the prime p, the dynamic system 
can be expressed in terms of the initial conditions. 


( 2 ) 


with 


n, 


A k = diag 


XPk*) = M\p k • • • M\ pi ■ n g (p 0 # ) 
= R ■ A k ■ L ■ n g (p 0 *) 


n (.Pi - 2 )> n<»- 3 )’■ • ■ ’fife - j ~ i ) 
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Dividing by p — 2 at each iteration, we obtain the normalized population 
model for gaps 

(3) Wg(Pk*) = M\ p fc ■ ■ ■ M\ P1 ■ Wg(p 0 *) 

= R ■ A k ■ L ■ w g (po # ) 

and here 


A k 


diag 


i TT — ~ ^ TT — ~ ^ 
’ | Pi — 2 : lift-2 


tt Pi ~ J — 1 
’ 11 Pi - 2 


The matrix L is an upper triangular Pascal matrix, whose entries are 
the binomial coefficients, and R is an upper triangular matrix whose entries 
have alternating signs. 

( (-1 ) i+i (Ci) ^ i<3 \ (ili 1 ) if *<J 

Rij = < and Lij = < 

[ 0 if i > j [ 0 if i > j 


This dynamic system encompasses constellations as driving terms for 
gaps. To adapt the system to focus directly on the population of a constel¬ 
lation and its driving terms, we only have to shift the eigenvalues. 


2.2. Extending the model to constellations. The model for gaps relies 
on keeping track of the internal closures and external closures for the driving 
terms for gaps. For constellations, we need a richer concept than external 
closures, and so we define the boundary closures for a driving term for a 
constellation. 

Let s be a constellation of length j\ , s = gi, g2, ■ ■ ■, gj 1 ■ A driving term 
s for s will have the form s = s~i s ~2 ... in which |sj| = gi for each i. 

We call the closures within an Si the interior closures for s, and the 
exterior closures for each Sj are the boundary closures for s. Interior closures 
preserve the copy of s as a driving term for s but of shorter length. The 
j i + 1 boundary closures remove the copy of s from being a driving term for 

s. 


We are now ready to prove Theorem ll.il 

Proof, of Theorem M.il - Let s be a constellation of length j\ and of sum s|. 
Let po be a prime such that |s| < 2p\. 

Let n s j{p be the population of driving terms of length j for s in <?(P # ), 
with ji < j < J and p > po- 

Consider the recursion from Q{pk-T^) to G(pk#), with k > 1. Let s be 
a particular occurrence of a driving term of length j in G(pk- w). Since 
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|-S* | < 2 pk, for the pf ; copies of s initially created during step R2 of the 
recursion, the j + 1 closures of step R3 all occur in different copies. 

Of the pk initial copies of s, the boundary closures eliminate j i +1 copies 
as driving terms for s. The j — j\ interior closures produce j — j\ driving 
terms of length j — 1, and pk — j — 1 copies of s survive intact as driving 
terms of length j. 


We can express this in the dynamic system: 

n s (pk*) = M ji:J (p k ) ■ n s (pk-i*) 

Pk~ ji- 1 1 0 • 

0 Pk ~ h ~ 2 2 • 


0 

0 

J-jl 
Pk - J -1 


Note that Mj 1 -j(p) is a (J — j\ + 1) x (J — j\ + 1) matrix. 


n g (Pk- 




The matrix Mj v j(p ) has the same eigenvectors as the matrix for gaps (of 
size J — j i + 1). The eigenstructure for Mj i: j(p ) is given by 

Mj i:J [p) = R ■ A j V .j(p) ■ L 

in which R is still the upper triangular Pascal matrix with alternating signs, 
and L is the upper triangular Pascal matrix. These do not depend on the 
prime p. The eigenvalues are 

A j v j(p) = diag \p — ji — 1, p — j! — 2,. . n p - J - 1]. 


Since the eigenvectors do not depend on the prime p , the dynamic system 
can easily be expressed in terms of the initial conditions. 


(4) 


with 


A j 1 .j = diag 


n s {Pk*) = Mjr-J\ Pk ■ ■ ■ M \ P1 ■ n S (Po*) 
= R ■ A jv.J ■ L ■ n s(P0*) 


YliPi-h- 1 ), Ikp.-n-D,..., YliPi-J- 1 ) 


□ 


2.3. Normalizing the populations. In the large, the population of a con¬ 
stellation of length j i grows primarily by a factor of (p — j\ — 1). So all 
constellations of length j ultimately become more numerous than any con¬ 
stellation of length j + 1, and comparing the asymptotic populations of 
constellations of different lengths is thereby trivial. 

On the other hand, to determine the relative occurrence among constel¬ 
lations of a given length j i, we divide by the factor ( p — j\ — 1). We define 
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the normalized popidation of a constellation s of length j\ and driving terms 
up to length J as 


Wsip*) = J] 


k<?>ii+l 


Q ~ Ji - 1 


n s (p # ) = 


<t>h+i(P*) 


n s {p 




We introduce the functions 4>i(p#) = \\ P q>i {q — i) 


Definition. Let Q = (pq-i • • • q m be a product of distinct primes, with 
qi < q 2 < ■ ■ ■ < q m - We define 

MQ) = (<!] -*)• 

qj >i 

Note that (pi = (p, the Euler totient function, over the defined domain - 
products of distinct primes. For the primorials we have 

v 

Mp*)= n(*-o. 

q>i 


For the normalized populations, the dynamic system becomes 


(5) w s (p k # ) = R ■ K k h .j ■ L ■ w s (p 0 *) 

with 



diag 1, 


][• 

i Pi~ J l - 1 


n, • q 

ir 

V Pi- 3 1 - 1 


n 


pi - J -1 

Pi - ji -1 


For these normalized populations, the dominant eigenvalue is now 1. 


For gaps, this normalization corresponds nicely to taking the ratio of the 
population of the gap under consideration to the population of the gap g = 2. 
For j i = 2, we can again interpret this normalization as a ratio, in this case 
to the populations of constellations s = 24 or s = 42. These constellations 
have no driving terms longer than ji, and U 2 a{p^) = ^ 42 {p^), so there 
is no ambiguity. This correlation between the normalization and ratios to 
constellations consisting of 2’s and 4’s begins to break down at j\ = 3 and 
collapses completely at j\ = 6. For j\ = 3, we have two constellations to 
choose from, s = 242 or s = 424, and there is ambiguity because 74424 (p^) = 
2 • U 242 (p^)- For j\ = 4 we have s = 2424 and s = 4242, and for j\ = 5 
we have s = 42424. However, by looking at we see that there are no 

constellations consisting only of 2’s and 4’s for length j\ > 6. 


The normalization does not necessarily provide ratios of the populations 
of constellations of length j\ to the population of a known constellation. 
Instead, it provides relative populations to a (perhaps hypothetical) con¬ 
stellation without driving terms other than the constellation itself and with 
population cf>j 1+ i(p #). 
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For example, for j\ = 3 we have po = 5 and po — ji — 1 = 1- The 
symmetric constellation s = 242 has w^ 2 = 1; it has no additional driving 
terms, it is symmetric, but it occurs at the middle of the cycle £?(5^). In 
contrast, the symmetric constellation s = 424 has rFm = 2, even though it 
too has no additional driving terms. 


2.4. Relative occurrence in the large. Equation (JSJ) applies under the 
same conditions as Theorem o - that s is a constellation of length j\ and 
driving terms up to length J, and po is a prime such that | s \ < 2p\. 

If we can get a count n s (po #) for s and all of its driving terms in G(po # ) 
with |s| < 2pi, then we can apply the eigenstructure to obtain the asymp¬ 
totic number of occurrences of s relative to other constellations of length 
ji as p k —5- oo: 

wf = lim w sjl (pk*) 

p k -KX> 

V n s,ji{Pk 

PK <j> jl + l( Pk #) 

= L 1 -w s {pq*) 

For this asymptotic result, we can soften the requirement of needing to 
work with po such that |s| < 2pi. For the full dynamic system to apply, this 
condition guarantees that for each occurrence s of a driving term for s in 
Q{pk- 1 ^), the closures in forming G(p k #) occur in distinct copies of s. 

Since L\ = [1,1,..., 1], to calculate we only need to know that the 
j i + 1 closures that remove a copy of s from being a driving term occur in 
distinct copies. Other closures may occur inside those copies as well, but 
we need to know that under the recursion p — ji — 1 copies of s survive as 
driving terms for s. 

Theorem 2.1. Let s be a constellation of j\ gaps 

s = 9i, 92 , 9 n ■ 

Let po be the highest prime that divides any of the intervals gi + ■ ■ ■ + gj with 
1 < i < j < ji■ Then for any p > po, 

= L { ■ w s (p*). 


Proof. Let s be a driving term for s. From G(po#) on > the closures that 
remove copies of s from being driving terms for s all occur in distinct copies. 
Thus for any subsequent p^, 


n s (pk*) = (Pk ~ ji - 1) n s (p k -i*) 
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We cannot be certain about the lengths of all of the copies that survive as 
driving terms, but we do know how the total population grows - by exactly 
the factor p^ — j\ — 1. and so for all pk > Po, 

wf = Li- w s (p k *). 

□ 


3. Specific constellations 


We now apply the population model to a few specific constellations. Pur¬ 
suing the twin prime conjecture, we study how constellations that corre¬ 
spond to twin primes persist in Eratosthenes sieve. One of these, s = 242, 
has been studied extensively before m mm-, this constellation corre¬ 
sponds to prime quadruplets, two pairs of twin primes separated by a gap 
of 4. 

The other family of constellations that we study here are those corre¬ 
sponding to consecutive primes in arithmetic progression. These constella¬ 
tions consist of the same gap repeated, such as s = 66 or s = 666. 


3.1. Constellations related to twin primes. Twin primes correspond 
to a gap g = 2. Here we are studying how gaps g = 2 arise in Eratosthenes 
sieve, and we do not address how many of these might survive the sieve to 
become gaps between primes (in this case between twin primes). 

The gap g = 2 arises in several interesting constellations. The first, 
s = 242, corresponds to prime quadruplets. Prime quadruplets are the 
densest occurrence of four primes in the large, two pairs of twin primes 
separated by a gap of 4. The constellation s = 242 has no additional driving 
terms, j\ = J = 3. We could use po = 3, for which n242(3^) = [1]. We have 

W 2A2 = 1 - 


3.1.1. Constellation s = 2,10,2. The next constellation we consider is s = 
2 ,10,2, which corresponds to two pairs of twin primes separated by a gap 
of 10. Here j± = 3 and J = 4, for the driving terms 2642 and 2462. Using 
Po = 7, we count n.2,io,2(7^) = [2,6] for two occurrences of s = 2,10,2 and 
three occurrences each of the driving terms 2642 and 2462. 


W2,10,2(7 # ) 

and 


1 

<M7 # ) 



w 2^10,2 — 8/3 

This means that as p^ —> oo, the number of occurrences of s = 2,10, 2 
in the sieve approaches 8/3 times the number of occurrences of the con¬ 
stellation 242. Remember that these weights w°° are relative only to other 
constellations of the same length j\. 
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3.1.2. Constellation s = 2,10,2,10,2. The constellation s = 2,10,2,10,2 
corresponds to three pairs of twin primes with gaps of g = 10 separating 
them. This constellation also contains two overlapping copies of 2,10, 2, and 
two overlapping driving terms for the constellation 12,12 to which we will 
return when we look at arithmetic progressions. 


For s = 2,10, 2,10, 2, we have j\ = 5 and J 
initial conditions we have to use Q{ 13^). 


= 7. Since |s| = 26, for 


n 2 , 10,2,10,2(13^) — 

^2,10,2,10,2(13^) = 


52 

44 

48 

1 


0 6 (13#) 


' 52 ' 

1 

' 52 ' 

44 

44 

48 

35 

48 


and 

^ 2 ) 10 , 2 , 10,2 = 144/35. 


So among constellations of length 5, the constellation s = 2,10, 2,10, 2 oc¬ 
curs with a relative frequency of 144/35. For length j\ = 5, we can use the 
constellation 42424 as a reference. The population model shows that in the 
large, the constellation s = 2,10, 2,10, 2 occurs over four times as frequently 
as the constellation 42424. 


3.1.3. Constellation s = 2,10, 2,10, 2,4, 2,10, 2,10, 2. Along this line of in¬ 
quiry into constellations that contain several 2s, we observe that the follow¬ 
ing constellation occurs in t/(13^): 

s = 2,10,2,10,2,4,2,10,2,10,2. 

This is two copies of 2,10, 2,10, 2 separated by a gap of 4. This corresponds 
to six pairs of twin primes, or twelve primes total, occurring in an interval 
of |s| = 56. 

For this constellation s, we have j\ = 11 and J = 13. It would theo¬ 
retically be possible for each of the 10’s to be produced through closures, 
so there could be driving terms of length up to 15, but when we inspect 
Q( 13 # ), we find that there are two copies of s, ten driving terms of length 
12, twelve driving terms of length 13, and no driving terms of length 14 or 
15. 


Since |s| = 56, to use Theorem 11.11 we would need to use po = 23 to 
employ the full dynamic system. However, to obtain the asymptotic results 
of Theorem 0 we can use po = 13. We calculate wf 3 = 24. This means 
that relative to a constellation of length 11 with one occurrence in C/(13^) 
and no additional driving terms, as the sieve continues the constellation s 
will occur approximately 24 times as often. 
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s 

|s| 

ji 

J 

Po 

n(p 0 *) 

w°° 

242 

8 

3 

3 

5 

[1] 

1 

424 

10 

3 

3 

5 

[2] 

2 

2,10,2 

14 

3 

4 

7 

[2,6] 

8/3 

42424 

16 

5 

5 

7 

[1] 

1 

2,10,2,10,2 

26 

5 

7 

13 

[52,44,48] 

144/35 

2,10,2,10,2,4,2,10,2,10,2 

56 

11 

13 

13 

[2,10,12] 

24 

66 

12 

2 

4 

5 

[0,2,2] 

2 

12,12 

24 

2 

6 

11 

[0,2,20,48,58] 

2 

666 

18 

3 

5 

7 

[0,4,2] 

2 


Table 1. Table of initial conditions and parameters for a 
few representative constellations. The population of a con¬ 
stellation of length j i grows primarily by a factor of p—j± — 1. 


3.2. Consecutive primes in arithmetic progression. A sequence of 
j i +1 consecutive primes in arithmetic progression corresponds to a constel¬ 
lation of ji identical gaps g. By considering residues, we easily see that for 
a sequence of j\ + 1 primes in arithmetic progression, g must be divisible by 
every prime p < j\ + 1. So for three consecutive primes in arithmetic pro¬ 
gression, the minimal constellation is s = 66. For four consecutive primes in 
arithmetic progression, the minimal constellation is s = 666, and then for an 
arithmetic progression of five consecutive primes the minimal constellation 
is s = 30,30,30,30. 

3.2.1. Constellation s = 66. Let’s now calculate the population of the con¬ 
stellation s = 66. This constellation corresponds to three consecutive primes 
in arithmetic progression: p, p + 6, p + 12. Since |s| = 12, we can still use 
po = 5. In g{ 5 # ) = 64242462, we observe the following initial conditions for 
s = 66: 

' 0 ' 

« 66 ( 5 # ) = 2 . 

2 

For the first entry, we don’t yet have any occurrences of s = 66. We do have 
two driving terms of length three: 642 and 246; and two driving terms of 
length four: 4242 and 2424. 


^66 (5 # ) 

and 

^66 



2. 
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3.2.2. Constellation s = 12,12. The constellation s = 12,12 also corre¬ 
sponds to three consecutive primes in arithmetic progression: p, p + 12, p + 
24. Since |s| = 24, to apply the full dynamic system of Theorem 11.11 we 
use po = 11. In £?(11#), we calculate the following, for driving terms from 
length ji = 2 to J = 6: 



0 


0 


2 


2 

ni2,i 2 (H # ) = 

20 

and Wi 2 ,i 2 (ll # ) - 

20 


48 


48 


58 


58 


So for s = 12,12 we have 


^oo 

Wl2 ’ 12 “ 8-4-2 


= 2 . 


It is interesting that the asymptotic relative population of s = 12,12 is the 
same as for the constellation s = 66. 


3.2.3. Constellation s = 666. The constellation s = 666 is the smallest 
constellation corresponding to four consecutive primes in arithmetic pro¬ 
gression. Since |s| = 18, we can use po = 7- In G(T&), we have the following 
initial conditions for s = 666, for driving terms from length j\ = 3 to J = 5: 


^666 ( 7 # ) 


" 0 " 

A 


" 0 " 

A 

2 

ana j 

2 


For s = 666 we have 

w 666 = 2 = 2 - 

It is interesting that we again have the asymptotic relative population of 
wf° = 2, although here it is relative to constellations of length j\ = 3. 


4. POLIGNAC RESULT FOR ARITHMETIC PROGRESSIONS 


Here we establish a strong Polignac result on arithmetic progressions. In 
[8j we established the following analogue of the Polignac conjecture: 


Theorem 4.1. Let g = 2n and let Q = q\ ■ ■ • q rn be the product of the 
distinct prime factors of g. Then g arises as a gap in Eratosthenes sieve 
and asymptotic ratio 


w 


OO 

9 


n 

2<! Qi 


qi - l 
qi- 2 ' 


For a gap g = 2n, we can now show that not only does the gap occur in 
Eratosthenes sieve, but that every feasible repetition of g as a constellation 
s = g,... , g occurs in the sieve. We make this precise below. 
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Above, we calculated the occurrences for a few small examples. We 
cannot perform the brute force calculation for five primes in arithmetic 
progression. The minimal constellation we would need to consider is 
s = 30, 30, 30, 30. To apply Theorem 11.11 we would have to use po = 57. 
However, we could obtain asymptotic results via Theorem 12.11 with po = 5. 
The steps we would take in order to apply Theorem 12.11 to this constellation 
can be generalized to prove Theorem 14.21 below. 

Definition. Let g be a gap and let pk be the largest prime such that 
Pk#\g- Let s = g, ..., g be a repetition of g of length j\. Then the constel¬ 
lation s is feasible iff j\ < Pk+i — 1. 

If it survives subsequent stages of the sieve, a repetition s of length j\ 
corresponds to j\ + 1 consecutive primes in arithmetic progression. 

Theorem 4.2. Let g be an even number, and let Q = q± - ■ ■ q m be the 
product of the distinct prime factors of g, with q\ < ... < q m . Let s be a 
feasible repetition of g of length j\. Then s occurs in Eratosthenes sieve with 
asymptotic weiqht 

wOQ MQ) 

9 ’-’ 9 tn+1 (QY 

Proof. We start in the cycle of gaps G{Q), in which Q may not be a primorial. 
G(Q) consists of <f>(Q) gaps that sum to Q. So we can start with any of the 
4>(Q) gaps and continue through the cycle g/Q times, and this concatenation 
is a driving term for g. 

We repeat this concatenation of cycles G{Q), to identify the driving terms 
for s = g ,..., g, a feasible repetition of the gap g of length j\. Starting with 
any gap in G{Q), we continue through G(Q) for j\ ■ g/Q complete cycles, 
and this is a driving term s for s. Since we can start from any gap in G(Q), 
we have 4>(Q) driving terms for s. 

Now that we have seeded the construction, we move from G(Q) back into 
the cycles of gaps for the primorials, G{q m#)- 

Case 1. Q is itself a primorial. In this case Q = and the total 

number of driving terms for s equals L\ ■ n s (q m = <f>(Q ). From this, we 
calculate the asymptotic ratio 

oo = ^(gm # ) = MQ) 

Ws (Ph+liQm*) <An+i(<2)’ 

Case 2. Q is not a primorial. Let p\ < p 2 < ■ ■ ■ < Pk be the primes 
less than q m that are not factors of Q. Let Qo = Q. For i = 1 ,..., k, let 

Qi — Qi—1 ' Pi- 

Let s be a driving term for s in G(Qi-i). Under the recursion to create 
G(Qi), we create pi copies of s in step R2. Then in step R3 we close gaps 
as indicated by the element wise product pi * G(Qi- 1 ). 
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The driving term s is composed of j\ driving terms for g. 

S <$2 • • • ^j\ 

in which each |sj| = g. 

The ji + 1 boundary closures will eliminate copies of s from being a 
driving term for s. All of the intervals in Theorem 12.11 are multiples of g, 
and these multiples themselves have prime factors entirely divisible by the 
factors of g. That is, if p& is the largest primorial that divides g , then for 
all 2 < j < ji, by the feasibility of s, the prime factors of j are factors of 

p*. 

Since pi J(g, all of the boundary closures for s occur in different copies of 
s. Thus, of the pi initial copies of s, j\ + 1 are removed as driving terms for 

s, and the other pi — j\ — 1 copies remain as driving terms of some length 

but no longer than the length of s. 

We don’t know the lengths of the driving terms for s in G(Qi), but we 
do know the total population: 

• n s (Qj) — (pj j i 1) ■ L\ • n s (Qi—i) 

= (Pi - Ji - l)(Pi-i — Ji — 1) — (pi - ji - 1) • Li ■ n s {Q 0 ) 

= (Pi ~ ji ~ l)(Pi-i - ji ~ !) • •' (Pi - ji ~ 1) • 4>{Q)- 

Continuing this construction, we have Qk = from which 

L\ ■ n s (Q k ) = ( pk - j i - l)(pfc-i - ji ~ 1) ■ ■ ■ (pi - Ji - 1) • &{Q), 
and the asymptotic ratio is 

? ,,o° _ \ _. T ■ n (a #) = 

s 1 s[qm } <t>ji+l(QY 

□ 

Thus every feasible repetition of a gap g occurs in Eratosthenes sieve. 
Compared to other constellations of length j]. the repetition of g of length 
j i has asymptotic weight (f>i(Q)/(f>j 1+ i(Q). To the extent that these con¬ 
stellations survive the sieve, they correspond to ji + 1 consecutive primes in 
arithmetic progression. 


5. Conclusion 

We have adapted the dynamic model for populations of gaps, as devel¬ 
oped in [8j, to apply to populations of constellations. These models for 
constellations all have the same left and right eigenvectors. The matrix of 
right eigenvectors R is the upper triangular Pascal matrix with alternating 
signs, and the matrix of left eigenvectors L is the upper triangular Pascal 
matrix. The eigenvalues depend on the prime and on the length of the 
constellation, but these otherwise have the same form. 
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We have calculated the populations for certain specific constellations con¬ 
taining the gap 2. These correspond to pairs of twin primes in various 
constellations. 

We also calculated the populations for a few constellations that are the 
repetition of a gap g. These repetitions correspond to consecutive primes 
in arithmetic progression. We are able to establish a strong Polignac result 
for arithmetic progressions: for any g = 2 n, every feasible repetition of g 
occurs as a constellation in Eratosthenes sieve. The asymptotic ratio for 
the occurrence of the repetition of g of length j\ depends only on the prime 
factors of g and the length j\. 
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